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ABSTRACT 

We present a statistical parallax analysis of low-mass dwarfs from the Sloan Digital Sky Survey 
(SDSS). We calculate absolute r-band magnitudes (M r ) as a function of color and spectral type, and 
investigate changes in M r with location in the Milky Way. We find that magnetically active M dwarfs 
are intrinsically brighter in M r than their inactive counterparts at the same color or spectral type. 
Metallicity, as traced by the proxy £, also affects M r , with metal poor stars having fainter absolute 
magnitudes than higher metallicity M dwarfs at the same color or spectral type. Additionally, we 
measure the velocity ellipsoid and solar reflex motion for each subsample of M dwarfs. We find good 
agreement between our measured solar peculiar motion and previous results for similar populations, as 
well as some evidence for differing motions of early and late M type populations in U and W velocities 
that cannot be attributed to asymmetric drift. The reflex solar motion and the velocity dispersions 
both show that younger populations, as traced by magnetic activity and location near the Galactic 
plane, have experienced less dynamical heating. We introduce a new parameter, the independent 
position altitude (IPA) , to investigate populations as a function of vertical height from the Galactic 
plane. M dwarfs at all types exhibit an increase in velocity dispersion when analyzed in comparable 
IPA subgroups. 



1. INTRODUCTION 

Low-mass dwarfs (0.08 M G < M < 0.8 Mq) 
dominate the stellar content of the Milky Way 
([Bochanski et al.l 12010ft . However, the study of these 
stars in large numbers was only recently realized with 
the advent of surveys such as the Sloan Digital Sky 
Survey (SDSS; lYork et all 12000ft and the Two Micron 
All Sky Survey (2MASS: ISkrutskie et all 12006ft . These 
surveys, with their precise, multi-band photometry and 
(in the case of SDSS) spectroscopic coverage have led 
to observational catalogs of unprecedented size. Cur- 
rently, the larges t spectroscopic database of M dwarfs 
()West et al.ll2011l hereafter Paper I) has produced spec- 
tral types, radial velocities (RVs) and chromospheric ac- 
tivity estimates (as traced by Balmer series emission) 
for over 70,000 M dwarfs. SDSS photometry of mil- 
lions of stars has been used to measur e the field luminos- 
ity and mass function s of M dwarfs (jCovev et all [2008: 
iBochanski et al.ll2010ft the structure of the Milky Way's 
thin and thick disks (|Juric et al.l 120081 : IBochanski et al.l 
12010ft . and th e properties of wid e, common proper mo- 
tion binaries (jDhitaretaD[20i3). 

The ubiquity and long main s e quenc e lifetimes 
of low-mass stars (jLaughlin et al.1 11997ft make M 
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dwarfs ideal trace rs of nearby Galactic st r ucture and 
kinem a tics (e.g.. iReid fc Maiewskil 119931: IReid et all 
1991 IBochanski et all I2007at iFuchs et al.l 120091: 
Bochanski et al. 120101 ) . Despite the recent advances 
enabling the study of these stars in a broad Galactic 
context, one fundamental parameter remains elusive: 
distance. Accurate trigonometric parallaxes, and the 
resulting absolute magnitude and distance measure- 
ments, remain difficult to obtain due to the intrinsic 
faintness (L < 0.05L©) of low-mass stars. The largest 
trigonometric parallax s urvey to date was performed b y 
the Hipparcos satellite ()ESAI 119971: Ivan Leeuwenl 12007ft . 
However, due to the relatively bright magnitude limit of 
the Hipparcos sample (V < 12), M dwarfs observed by 
Hipparcos are saturated in SDSS photometry, hindering 
the construction of a color-absolute magnitude relation 
(CMR) in native SDSS filters. A few parall ax estimates 
of SPSS M dwarfs hav e been obtained (e.g.. lDahn et all 
l2002trVrba et alll2004ft . but these measurements are ob- 
servational!/ taxing and take years to comp lete, and are 
thus limited in number (<~ 20). Current (jRiedel et all 
120101 : Faherty, private comm unication) and future 
parallax studies, such as GAIA (|Perrvman et all [2001ft 
will add valuable observations to this regime. 

Since trigonometric parallaxes for the vast majority of 
SDSS M dwarfs are not available, alternate means to 
estimate thei r abso lute magnitudes must be employed. 
IHawlev et all (|200l and I West et all ((2005ft used SDSS- 
2MASS cross-matching to create Mj-spectral type re- 
lations, and combined with average colors as a function 
of spectral type, to construct SDSS C MRs. Golimowsk i 
et al. (private communication) and iBochanslril (2008) 
used u'g'r'i'z' photometry of ~ 200 nearby M dwarfs 
with accurate trigonometric parallaxes, and transformed 
from the u'g'r'i'z' system of the SDSS photometric tele- 
scope t o the native SPSS ugri z system using the rela- 
tions of iDavenport et al.l (|2007ft . Other CMRs applica- 
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ble for warmer stars i n the SPSS filter- set have been 
derived using bin aries (iSesar et al. [2008) and averages 
of existing CMRs (jJuric et al.l l2008). Using a kinematic 
model for motions of stars in the Milky Way, iBond et al.1 
(2010) derived a distance scale and CMR for main se- 
quence stars observed by SDSS. Other studies have used 
2MASS-SDSS color transf ormations (iBilir et all I2009T) 
and synthetic photometry (|Covev et al.ll2007t ) to derive 
CMRs in the ugriz filter system. 

One method for constructing CMRs in the native 
SDSS system that has not yet been utilized is the clas- 
sical statistical parallax method. Using proper mo- 
tions and radial velocities of a large number stars 
from a homogeneous population, an estimate of the av- 
erage absolute magnitude, the velocity ellipsoid, and 
the Sun's peculiar motion can be obtained. This 
method has a rich histor y in the astron omical litera- 
ture (e^J^nHerSJTlel OubeXjonei [Till IMurravl 
119831: IPopowski fe Gouldl [19981), and has previously 
been employed for RR Lyraes (e.g., lHawlev et all 119861: 
Strugnell etaU 119861: lLavden et all 119961: iFernlev et all 
19981: IPopowski fc Gouldl 11998ft and Cepheids (e.g., 
Wilson et al.lll991[ ). 



In addition to absolute magnitude estimates, the sta- 
tistical parallax technique provides a measurement of the 
velocity ellipsoid of the stellar population and the reflex 
motion of the Sun. The velocity ellipsoid is described 
by dispersions and directions along three principal axes. 
Low-mass sta rs make excellent t racers of the local Galac- 
tic potential (iReid et al.1 119971 : iBochanski et al.1 12007a! : 
iFuchs et al.l 120091) as constrained by the observed veloc- 
ity dispersions. The solar peculiar motion determina- 
tion complements studies usi ng Hipparcos proper mo - 
tions of nearby, young stars dDehnen fc Binnevl I l998). 
SDSS observations of M dwarfs (|Fuchs et al.l 120091 ) and 
the Palomar- Michi gan State University survey of nearb y 
M dwarfs (PMSU: iReid et al.lll995l: lHawlev et al.l[l996h . 

In this paper, we use the large spectroscopic sample 
of M dwarfs described in iPaper II to determine absolute 
magnitudes in the native SDSS filters and examine kine- 
matics using the statistical parallax method. The obser- 
vations employed in our study are described in Sj2] We 
detail our statistical parallax analysis in SJ3] Our results 
are presented in 21 with conclusions following in §5] 

2. OBSERVATIONS 

Accurate and precise photometry, astrometry and ve- 
locities are required for statistical parallax analysis. The 
data w ere obtained from the latest SDSS data release 
(DR7; lAbazaiian et al.l 120091) which contains photom- 
etry over nearly 10,000 sq. deg. down to faint magn i- 
tudes (r ~22) in five filters (ugriz, iFukugita et alf l996). 
When sky conditions at Apache Point Observatory were 
not photometric, the SDSS operated in a spectroscopic 
mode. SDSS photometry was used to tar get objects for 
spect roscopic follo wup, primarily galax ies (|Strauss et al.l 
2002) and quasars (|Richards et al]|20"02j ). Approximately 
460,000 stellar spectra were also obtained both as tar- 
geted and serendipitous observations. Twin fiber-fed 
spectrographs collected 640 observations simultaneously, 
with individual 15-20-minute exposures being co-added 
for a typical total exposure time of ~ 45 minutes. These 
medium-resolution (R ~ 2,000) spectra cover the entire 
optical bandpass (3800 - 9200 A). 



The absolute ast rometric precisio n of SDSS is < 0.1" 
in each coordinate (|Pier et al.ll2003D . Proper motions for 
SDSS obj ects were calculated by m atching to the USNO- 
B survey (jMunn et alJl2004l 120081 ). The proper motions 
have a baseline of ~ 50 years, and a typical precision of 
~ 3 mas yr - 1 in right ascension and declination. 

In IPaper 1 we described the SDSS DR7 spectra that 
we compiled into the largest spectroscopic catalog of 
M dwarfs ever constructed, containing ~ 70,000 stars. 
Briefly, the catalog was obtained by color selection of 
stars with r — i > 0.42 and i — z > 0.24, befor e 
correcting for Galactic reddening (jSchlcgel ct al. 1998). 
Each star was processe d with the Hammer IDL pack- 
age (jCovev et al.l 12007ft . The Hammer provides spec- 
tral type estimates and measures a number of spec- 
tral features, including Ha equivalent width and vari- 
ous TiO and CaH bandhead strengths. Each spectrum 
was then visually inspected and non-M dwarf contam- 
inants were culled from the catalog. Spectral types 
were determined by com paring to the SDSS templates 
(|Bochanski et al.ll2007bD . The RVs were measured by 
cross-correlati ng; each spectrum with the appropriate 
template from Bochanski et al. (2007bj), with a precision 
of ~ 7 km s -1 . See IPaper 11 for more detailed informa- 
tion about the selection and characteristics of the DR7 
catalog. 

For the sample used in the statistical parallax analysis, 
we required that stars have well-measured proper mo- 
tions and radial veloci ties, and precise photometry (see 
Bochansk i et al.l (12010ft fo r infor mation on photometric 
flag cuts, an d IDhital et afl <l27TT?^ for our kinematic qual- 
ity flags). Our quality cuts reduced the original sample 
to 40,963 stars. 

3. METHOD 

We employed the maximum likelihood formulation 
of classical statis tical parallax analysis as presented 
by IMurravl (119831) and used pre v iously to analyze RR 
Lyrae stars (lHawlev et al.l 119861: iStrugnell et~aLl 119861: 
19961 |Fernky_eLaJL| UMM) and Cepheids 
1991). A variation of the model was used 
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to investigate the k inematics of the nea rby M dwarfs from 
the PMSU sample ([Hawlev et al.lll996ft . Briefly, this sta- 
tistical parallax method models the velocity distribution 
of a homogeneous stellar population with nine kinematic 
parameters, including the three velocities of the reflex 
solar motion, and the three directions and three disper- 
sions of the velocity ellipsoid, which describes the ran- 
dom and peculiar velocities of the population. The stars 
are assumed to have an absolute magnitude with some 
intrinsic dispersion, and these two additional parameters 
determine the distance used to transform the observed 
proper motions into transverse velocities. The eleven pa- 
rameters used in the model are described in Table [TJ 

The observational data needed for the analysis are the 
position, apparent magnitude (corrected for Galactic ex- 
tinction), proper motions and radial velocity for each star 
in the sample. The uncertainties in the data are used to 
assign appropriate weights in the solution. We solved for 
the eleven parameters in the model simultaneously by 
maximizi ng the likelihood using geometric sim plex opti- 
mization (|Nelder fc Meadl H9651: IDanielsl li~978ft . Uncer- 
tainties in the parameters were determined by numerical 
computation of the derivative of each parameter indi- 
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TABLE 1 

Statistical Parallax Fit Parameters 



Parameter 


Units 


Description 


<W 


km s — 1 


Velocity dispersion in radial direction 


a v 


km s — 1 


Velocity dispersion in orbital direction 


cr w 


km s — 1 


Velocity dispersion in vertical direction 


r, (f>, z 


radians 


Orientation of velocity ellipsoid 


U 


km s~ 1 


Solar peculiar motion (r) 


V 


km s~ 1 


Solar peculiar motion (</>) 


w 


km s — 1 


Solar peculiar motion (2) 


k 




Distance scale 






Dispersion in distance scale 



vidually, while keeping all other parameters fixed. The 
maximum likelihood equations and simpl ex method are 
described in detail in ' Hawlev et al.l ()1986[ ). 

With any multi-parameter fit, sensible constraints are 
necessary to ensure the model does not prefer an unreal- 
istic portion of \ 2 space. For this analysis, we fixed one 
parameter, o~k- This parameter is related to the spread in 
absolute magnitude , <jm by the following equation from 
((19861) : 
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cV = log 10 [l + (7 fe 2 /(l + fc) 2 ]/(0.04 In 10) (1) 

where o~m is the spread in absolute magnitude for a given 
luminosity bin, and <7fc is the spread in k which is a dis- 
tance scale parameter. The absolute magnitude is related 
to k through the following equation: 



Mr 



51og 10 (l + fc) + M j4 -0.1 lnl0(<7 M ) 



(2) 



where Ma is an initial estimate of the absolute magni- 
tude. We fixed Uk at four values: 0.05, 0.1, 0.2 and 
0.3. Usually k ~ 0.01 — 0.02, so at — 0.2 corresponds 
to <tm ~ 0.4 which is the typi cal dispersion in absol ute 
magnitude for low-mass stars (jBochanski et alJ l2010). 



3.1. Constructing Subsamples for Analysis 

Our spectroscopic sample is much larger than was 
available for previous statistical parallax studies, which 
contained observations of a few hundred or less objects. 
Thus, we were able to divide the sample into smaller sub- 
samples, selected on color, spectral type, magnetic activ- 
ity (a s traced by Ha), m etallicity (using the £ parame- 
ter of iLepine et al.ll2007l ). position on the sky (northern 
and southern Galactic hemispheres) , dis tance (estimated 
from the M r , r — z CMR relation from iBochanski et al.l 
2010), and Galactic height (discussed further below). We 
required at least 100 stars in each subsample to ensure 
a meaningful solution. Table [2] lists the properties of the 
various subsamples. 

The color, spectral type, metallicity and magnetic ac- 
tivity subsamples were chosen to explore intrinsic varia- 
tions within the M dwarf population, as traced by their 
absolute magnitudes and kinematics. The position and 
distance subsets were examined to ensure that we were 
not biased by extinction or proper motion limits. While 
the majority of our sample is contained at high Galactic 
latitudes (b ~ 50°) in low extinction regions (the median 
extinction in r is 0.07 mags), some sightlines to the south 
(through the Plane) may be affected by extinction and 



reddening. Current efforts to model the 3-D dust distri- 
bution within the SDSS footprint (Jones et al., private 
communication) will improve the extinction and redden- 
ing estimates for these stars. 

Our prope r motion catalog, described in iMunn et al.1 
(|2004L |2008| ) has a precision limit of 3 mas yr -1 . For 
a typical thin disk star with a transverse velocity of 
~ 10 km s _1 , this limit corresponds to a limiting dis- 
tance of ~ 700 pc. Beyond this distance, only stars 
with larger transverse motions would be measured with 
precise proper motions, artificially biasing the measured 
velocity dispersion. Thus, two distance limits were ex- 
plored: stars with d < 500 pc and d < 700 pc. The 
500 pc cut corresponds to a more conservative estimate 
of the proper motion precision. The effects of these dis- 
tance limits are shown in Figure[TJ Early-type M dwarfs, 
which are found at larger distances, are biased towards 
larger velocity dispersions if a distance limit is not en- 
forced. The 500 pc and 700 pc limits give similar results, 
so we adopt a d < 700 pc limit for the remainder of this 
study. Further kinematic analysis is discussed in £14.21 
The d < 700 pc cut limits the total number of stars in 
our sample to 22,542. In Table [3] we report the median, 
minimum and maximum heights above the Plane for the 
stars in our d < 700 pc sample. The heights were mea- 
sured usin g distances from the M r , r — z color-magnitude 
relation of Bocha nski et al.l (|2010[ ) . 
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Fig. 1. — UVW dispersions vs. spectral type for d < 500 pc 
(open triangles), d < 700 pc (open squares) and no distance limit 
(filled circles). Including more distant stars beyond our proper 
motion precision limit artificially inflates the measured dispersions 
for early— type stars. 



Finally, to probe changes in the kinematics and ab- 
solute magnitude as a function of height in the Milky 
Way disk, we divided each spectral type and color bin 
according to the independent position altitude (IPA) . At 
a given spectral type, bins in apparent magnitude will 
correspond to spherical shells in space. The IPA cuts are 
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slices in Galactic height, defined as: 

IPA = 10 m -- /5 sin& (3) 

where m r is the observed r-band magnitude and b is 
the Galactic latitude of the star. We stress that the 
IPA can only be used for relative comparisons within 
the same intrinsic luminosity bins. That is, MO and M6 
stars with the same IPA value will be located at very dif- 
ferent Galactic heights. The IPA is useful as it avoids the 
assumption of an absolute magnitude, allowing an inde- 
pendent determination for each IPA subset. It is based 
entirely on observable properties, and is somewh at anal- 
ogou s to the often-used reduced proper motion ([Luvtenl 
1925). 



3.2. Computational Method 

Following the method of lHawlev et al.l (|1986D each sub- 
sample of M dwarfs was analyzed using the following pre- 
scription. Ten runs were calculated for each dataset. The 
initial absolute magni tude estimates w e re de rived from 
the M r ,r-z CMR of iBochanski et al.l (|2010f ). The ab- 
solute magnitude estimate was updated after each run, 
with the output of the previous run being used as input 
for the next. For a given run, the simplex optimization 
iterated 5,000 times. Thus, for each subsample, 50,000 
iterations were computed. This ensured that the sim- 
plex operator had sufficient freedom to explore parameter 
space and converged to the best (maximum likelihood) 
answer. Typically, convergence was obtained after 25,000 
iterations. 

4. RESULTS 

The results of our analysis are detailed below. The 
outputs from the code returned three major results: the 
absolute magnitude of a given stellar sample, the Sun's 
velocity with respect to that sample, and the velocity 
ellipsoid that describes the peculiar motions of the sam- 
ple. Each of these results are examined, and compared 
to previous investigations. In the discussion, we use the 
(Tfc = 0.2 results as our fiducial value, since it corresponds 
to <jm t — 0.4 which is observed for nearby stars with 
trigonometric parallaxes in the M r , r—z color-magnitude 
diagram. 

4.1. Absolute Magnitudes 

Traditionally, colors or spectral types have been em- 
ployed to estimate absolute magnitude when trigono- 
metric parallaxes were not available. These are of- 
ten referred to as photometric and spectroscopic paral- 
laxes, respectively. In Figure [3J we show the spectro- 
scopic (left panel) and photometric (right panel) paral- 
lax relations for M dwa rfs, compared to pr e vious stud- 
ies in the uqriz system (lHawlev et al.l [20021 : iWest et"aTI 
120051 : iKraus fc Hillenbrand! I2007t ). The open circles are 
nearby low-mass sta rs with accurate trigonometric par- 
allax es, tabulated in IBochan ski (2008). Using the Ham- 
mer ()Covev et al.l 12007ft . we manually assigned a spec- 
tral type to each nearby star, using archi val spectra 
from multiple sources: the PMSU study (|Reid et al.l 



119951: lHawlev et al.lll~996l ). SDSS, the local star sample 
of ICruz fc Reicfl ([20021 ) and the Dwarf Archives websiteB 
While many of the spectral types agreed with previous 
results, we did find some disagreements at the level of ±2 
spectral classes. Therefore, we include updated spectral 
types for these nearby stars in Table [4] The left panel of 
Figure [2] highlights the large dispersion in M r as a func- 
tion of spectral type. The typical spread is ~ 0.5 mag, 
increasing to ~ 0.8 magnitudes near type M4, nearly 
double the spread for a given r — z color bin. Thus, we 
strongly suggest that spectral type should not be used 
as a tracer of absolute magnitude for M dwarfs. The use 
of half-interger spectral types may decrease the observed 
scatter in Figure [2] However, half-integer spectral types 
would require higher resolution and higher S/N spec- 
tra, which is not practical for this type of large survey 
sample. Furthermore, half-interger spectral type stan- 
dards are not defined f or the entire M dwarf sequence 
(jKirkpatrick et aD ll991). Finally, half-integer types will 
likely span a l arger co lor range than the color bins we 
have used (see iPaper II for a discussion of spectral types 
and mean colors). 

The left panel of Figure [2] also displays disagreement 
between our statistical parallax results (solid black, red 
and blue lines and filled circles) and the nearby star sam- 
ple at later types. This is due to a systematic color dif- 
ference within a given spectral type bin between the two 
samples. That is, the nearby star sample is systemat- 
ically redder than the SDSS stars at a given spectral 
type, possibly due to the small number of nearby stars 
in these bins. Thus, the SDSS stars have an absolute 
magnitude (at a given spectral type) that is appropriate 
for the bluer, more luminous stars of that type. This ef- 
fect disappears when r — z color is substituted for spectral 
type (Figure [21 right panel). 

The right panel of Figure [2] shows our statistical 
parallax M r ,r — z relation, along with the results of 
previous studies (lHawlev et al.l l200l IWest et al.l 120051 : 
IBochanski et al.l 120101 ). Overplotted with open circles 
are the absolute magnitudes and colors from the nearby 
star sample listed in Table @] Note that the dispersion 
in absolute magnitude as a function of color is signifi- 
cantly smaller than that shown in left panel of Figure [2j 
indicating that the r — z color is a much better tracer 
of absolute magnitude. Furthermore, the discrepancy in 
mean color for a given spectral type between the nearby 
star sample and SDSS stars is not evident. We reiterate 
that colors, rather than spectral types, are preferred for 
absolute magnitude and hence distance estimation for M 
dwarfs. 

4.1.1. Absolute Magnitude Variations and Magnetic Activity 

Due to the large spread in M~ r as a function of spectral 
type, we focus our absolute magnitude analysis on our 
color subsamples (i.e., the right panel of Figure [2]). Our 
results are shown in the solid black, red and blue lines 
and filled circles in Figure [2] There are a few notable 
trends. First, the bifurcated main sequence at blue colors 
(earlier spectral types) indicates that active, early-type 
M dwarfs are intrinsically more luminous than their in- 
active counterparts. Recent observations of eclip sing bi- 
naries (|L6pez-Moralesll2007l : iMorales et"alll2010[ ) and ac- 

6 Available at http://www.dwarfarchives.org 
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TABLE 2 

SDSS SlTBSAMPLES 



Initial Cut 



Bin Size 



Further Criteria 



Comments 



Spectral Type 
Color 



1 type bins 
0.3 mag bins 



Activity Ho active stars defined in F^arjer_I 

Hemisphere N and S Galactic Hemispheres 

C 0.25 bins 

IPA 10 bins per spectral type 

Distance < 500 pc k, < 700 pc 



Note. — Bins with fewer than 100 stars were not used in this analysis. 



TABLE 3 

Galactic Heights of SDSS M Dwarfs 





Median Height 


Min. Height 


Max Height 


Altars 


Spectral Type 


M0 


395 


20 


690 


1323 


Ml 


392 


22 


689 


1877 


M2 


337 


11 


693 


4700 


M3 


274 


14 


712 


6156 


M4 


214 


9 


694 


4784 


M5 


137 


9 


694 


1616 


M6 


108 


1 


423 


1333 


M7 


82 





230 


753 


r — z 


0.75 


402 


20 


690 


348 


1.05 


395 


23 


688 


1647 


1.35 


375 


16 


693 


2914 


1.65 


309 


11 


712 


5937 


1.95 


254 


13 


694 


5639 


2.25 


181 


9 


634 


2922 


2.55 


123 


11 


387 


818 


2.85 


118 


7 


304 


958 


3.15 


98 





221 


960 


3.45 


76 


2 


340 


296 


3.75 


58 


5 


126 


102 



Note. — All heights reported in pc and measured us ing dis- 
tances from the M r , r - z CMR of IBochanski et al.l (|2010T ). 



tive, single stars ()Berger et al.H2006l: IMorales et al.ll2008fl 
suggest that active stars possess larger radii than inac- 
tive dwarfs. However, possible changes in effective tem- 
perature and luminosity have no t been well constrain ed. 
Models pr edict both unchanged (IChabrier et al.ll2007h or 
decreased (jMullan fc MacDonaldll2001|) luminosity, while 
some observations suggest that active low-mass dwarfs 
are over-luminous compared to their mo del predictions 
(jMohantv et al.ll2009t iStelzer et aTll2010f ). A similar ef- 
fect was observed in the PMSU sample, where active M 
d warfs were brighter in V than inactive stars (Figure 4 
of lHawlev et al lll996[ ). The behavior evidenced in Fig- 
ure [2] suggests that active, early-type M dwarfs are ~ 1 
mag brighter in M r , which could be explained if they 
have larger radii than inactive stars at the same spec- 
tral type or color. Direct interferomctric measurements 
of stellar radii are necessary to independently quantify 
the amount the radius changes for an individual star, 
since our statistical analysis is performed on subsamples 
of stars that span a (small but nonzero) range of temper- 
ature. We note that many close binaries are unresolved 
in SDSS photometry. Binarity may increase activity in 
early M dwarfs in close orbits, and could explain some of 
the absolute magnitude differences for those bins. 



At blue colors (r — z < 2), the inactive and total sample 
absolute magnitudes fall below the mean locus of nearby 
trigonometric parallax stars. Since the nearby stars are 
composed of a mix of active and inactive stars, activity 
is likely not the only important effect. We therefore next 
investigate the effects of metallicity on M r and attempt 
to isolate them from those traced by chromospheric ac- 
tivity. 

4.1.2. Absolute Magnitude Variations and Metallicity 

While magnetic activity may affect the radius of a star, 
metallicity can alter the star's effective temperature and 
luminosity, manifesting as a change in absolute magni- 
tude. At a given color (or spectral type), stars with 
lower metallicities will exhibit fainter absolut e magni- 
tudes (ISandage fc Egged [19591 [An et al.l [2()08lf 7 l . With- 
out precise metallicities and trigonometric parallaxes for 
all the M dwarfs in our sample, we cannot directly exam- 
ine this effect. There are some promising metallicity es- 
timators for low-mass stars being developed in the near- 

7 Alternatively, at a given mass, stars with lower metallicities 
are bluer (hotter) and more luminous than their high metallically 
counterparts. 
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IR l| Johnson fc Appsll2009tlRoias-Avala et al.ll2010fl . but 
these are not applicable to SDSS spectroscopy. 

We therefore investigated the eff ects of metallicity us - 
ing the ( parameter, as defined bv iLepine et al.1 {2007). 
£ is a relative metallicity proxy which uses the relative 
strengths of CaH and TiO bands to discriminate between 
stars of different composition. It can be employed as a 
rough tracer of [Fe/H], with ( — 1 corresponding to solar 
metallicity and C = 0.4 corresponding to [Fe/H] ~ — 1 
(|Woolf et al.ll2009f ). We note that the [Fe/H]X relation 
is only calibrated over spectral types M0-M3 and suffers 
from large spreads near solar metallicity. In Figure [3j we 
plot the absolute magnitudes for two bins in £ as a func- 
tion of spectral type (left panel) and r — z color (right 
panel) for active (solid lines) and inactive (dashed lines) 
stars. Table [S] gives the M r results for the various r — z, 
activity and £ bins. The spread in each bin is a ~ 0.4 
mags, set by the choice of a% = 0.2. As with any mag - 
nitude limited survey (such as SDSS), |Ma!mgjfist| (1936) 
bias can arise. For surveys with complicated selection 
effects, a magnitude-inde pendent correction is preferred 
(jPopowski fc GouTdlll998f ). For an exponential Galactic 
disk a nd <jm, = 0-4, the M almquist correction is 0.15 
mags (Rcid & Hawlcy 2005). This correction has been 
applied to the values reported in Table [5J and is much 
smaller than the reported uncertainty. 

Figure [3] demonstrates that for stars with similar chro- 
mospheric properties, the lower metallicity ones have 
fainter absolute magnitudes, at the same color or spec- 
tral type. This is co nsistent with cl uster studies within 
the SDSS footprint (|An et al.l 12006) and trigonometric 
parallax studies of subdwarfs (|Reidlll997l ) . Active stars 
at the same metallicity are brighter than their inactive 
counterparts, but Figure [3] shows that both metallicity 
and activity are important for determining the luminos- 
ity of an individual star. At a fixed £, the difference in 
M r is consistent with Figure [2J with early-type active 
stars being ~ 1 mag brighter, and the difference dimin- 
ishing at later spectral types. A similar offset was also 
measured for the same stars in M g . The total SDSS sam- 
ple is overplotted in each figure with a solid black line. 
The early-type stars, which are seen at larger distances 
due to SDSS magnitude limits (see Tabled, fall near the 
lower metallicity, inactive loci, while the later type stars, 
seen closer to the Sun, are consistent with solar metallic- 
ities and active stars. This suggests that the low-mass 
dwarfs are tracing a metallicity gradient similar to the 
one observed in SP SS observations of higher-mass stars 
(jlvezic et a l. 2008) , and also explains why the SDSS sam- 
ple falls below the locus of nearby stars at early types 
(right panel of Figure [2]) . 

We also note that the activity-metallicity loci appear 
to converge near M5 (r — z ~ 2.8). This behavior is 
not well explained, but may be linked to the transition 
between a partially and fully convective stellar interior 
that occurs near that spectral type/color. Perhaps this 
transition, which alters the efficiency of energy transport 
in the star, also regulates the luminosity at the surface. 

Finally, we examined the IPA and hemisphere subsam- 
ples. No significant differences were observed in the spec- 
troscopic or photometric parallax relations with respect 
to position in the Galaxy. A trend with IPA was ex- 
pected, since this proxy for height probably traces metal- 
licity. The lack of a gradient is probably due to stars at 



various heights being scattered into the same IPA bin 
due to photometric errors. The agreement between the 
hemisphere subsamples indicates that the extinction cor- 
rection is adequate for the sample. 

TABLE 4 

Nearby Star Sample Spectral Types 



Name a a <5 a Sp. Type Source* 5 



GJ1002 


1.67972278 


-7.538550 


M6 


1 


GJ1025 


15.23465440 


-4.449231 


M5 


1 


G154.1 


18.12784543 


-16.998746 


M5 


1 


LHS1302 


27.76670103 


-6.117972 


M5 


1 


LHS1326 


30.56747119 


10.337209 


M6 


1 


LHS1358 


33.22773200 


0.004801 


M4 


1 


LHS1375 


34.12457669 


13.586802 


M6 


1 


GL109 


41.06474242 


25.523168 


M2 


1 


T832-10443 


43.10949326 


0.939522 


M9 


3 


LHS168 


48.34574221 


4.774824 


M5 


1 


GJ1065 


57.68434297 


-6.095079 


M3 


1 


G1169.1A 


67.79870688 


58.976686 


M4 


1 


LHS1723 


75.48913133 


-6.946128 


M4 


1 


G99-49 


90.01463326 


2.706515 


M4 


1 


LHS1809 


90.62144398 


49.865368 


M5 


1 


G1232 


96.17218963 


23.432802 


M4 


1 


G1251 


103.70407576 


33.268001 


M3 


1 


GJ1093 


104.87030750 


19.348539 


M5 


1 


BL Lyn 


112.98873721 


36.229615 


M4 


1 


2M0746+20 


116.67677307 


20.008842 


L0 


1 


GJ1105 


119.55303380 


41.303596 


M4 


1 


GJ2066 


124.03289059 


1.302435 


M2 


1 


GJ1111 


127.45563147 


26.775994 


M8 


1 


T213-2005 


155.36376953 


50.917938 


L0 


4 


LHS283 


158.86050669 


69.449381 


M4 


1 


G1445 


176.92332403 


78.691261 


M4 


1 


G1447 


176.93515470 


0.804025 


Ml 


1 


GJ1151 


177.73956422 


48.377048 


M4 


1 


GJ3693 


178.46963501 


6.998430 


M8 


2 


G1452.4 


178.73946569 


28.737442 


K7 


1 


G1455 


180.57503136 


28.587004 


M3 


1 


GJ1156 


184.74706541 


11.126184 


M6 


1 


G1463 


185.75015884 


64.030910 


M3 


1 


GJ1159A 


187.30915386 


53.545801 


M4 


1 


LHS2633 


191.75375224 


46.625855 


M2 


1 


G1493.1 


195.13932515 


5.685628 


M5 


1 


G1514 


202.49948063 


10.376790 


Ml 


1 


G1521 


204.85046387 


46.186687 


Ml 


1 


LHS2784 


205.68017011 


33.289742 


M4 


1 


LHS2884 


213.82055327 


45.014702 


M3 


1 


GJ3849 


217.17973328 


33.176899 


L0 


4 


G1552 


217.37332947 


15.533145 


M2 


1 


GJ3855 


217.65687658 


59.723624 


M8 


1 


G1555 


218.57030665 


-12.519293 


M4 


1 


2M1501+22 


225.28405762 


22.833836 


M9 


1 


LHS3018 


226.07657872 


60.384562 


Ml 


1 


G1581 


229.86143531 


-7.722272 


M3 


1 


G1585 


230.96299622 


17.465422 


M4 


1 


LHS3080 


232.97554057 


28.852619 


M4 


1 


G1609 


240.71170993 


20.588769 


M4 


1 


G1625 


246.35287094 


54.304172 


M2 


1 


G1628 


247.57513257 


-12.662992 


M4 


1 


LTT14949 


250.20368457 


36.316638 


M2 


1 


GL643 


253.85475375 


-8.322944 


M4 


1 


GJ1207 


254.27401671 


-4.349044 


M4 


1 


LHS3262 


255.84984701 


51.406478 


M4 


1 


GJ1209 


256.09302479 


16.931814 


M3 


1 


G1655 


256.78099478 


21.554112 


M3 


1 


LTT15087 


257.89507230 


38.442849 


M4 


1 


G1678.1 


262.59473243 


5.548422 


M0 


1 


G1686 


264.47272430 


18.592155 


Ml 


1 


G1694 


265.98321814 


43.378370 


M2 


1 


GJ1223 


270.69304642 


37.517084 


M5 


1 


G1701 


271.28155234 


-3.031352 


Ml 


1 


G11227 


275.61219859 


62.049964 


M5 


1 


G1720B 


278.86411519 


45.761550 


M4 


1 


G1725B 


280.69435149 


59.627614 


M3 


1 


G1729 


282.45607171 


-23.836125 


M4 


1 
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TABLE 4 — Continued 



Name « a <5 a Sp. Type Source b 



G1745A 


28fi 77322fl3fi 


20 887782 


M2 


x 


G1745B 


286 8(15(1671 6 


20 876758 


M2 


1 


G207-22 


288.12267278 


35.564413 


M2 


1 


G1752A 


289.22996750 


5.168470 


M3 


1 


GJ1235 


290.41104011 


20.867176 


M4 


1 


GJ1253 


306.52229530 


58.573236 


M5 


1 


G1809 


313.33253189 


62.153943 


Ml 


1 


LHS3713 


327.06371612 


27.927964 


M2 


1 


G1849 


332.41848642 


-4.640774 


M3 


1 


G1127-35 


337.19158157 


18.931698 


M0 


1 


G1867B 


339.68865820 


-20.614020 


M4 


1 


G1867AC 


339.68997030 


-20.620763 


M2 


1 


G1908 


357.30266119 


2.400989 


Ml 


1 



a and S are reported in decimal degrees and J2000 coordinates. 
b 1-PMSU IReid et all 199511 . 2-DwarfArchives.org, 3-SDSS llPaper IB . 
4- ICruz fc Reidl | [2002Tl 



TABLE 5 

Absolute Magnitudes of SDSS M Dwarfs 





c = 


1.1 


C = 


0.88 


r -z 


Active 


Inactive 


Active 


Inactive 


1.05 




8.8 






1.35 




9.1 




9.7 


1.65 


8.4 


9.8 


9.3 


10.4 


1.95 


9.0 


10.4 


10.3 


11.1 


2.25 


10.3 


11.3 


11.1 


12.0 


2.55 


11.2 


12.4 


12.5 




2.85 


13.2 


13.5 


13.6 




3.15 


13.9 


14.0 


13.9 




3.45 


14.9 









Note. — The typical uncertainty for the reported absolute magni- 
tudes is (TM r — 0.42, which is directly computed from EquationfTI 



4.2. Kinematics of M Dwarfs 

In addition to an absolute magnitude estimate, the 
analysis returns kinematic information for each subsam- 
ple. The interpretation of the kinematic results is com- 
plicated by the disjoint spatial volumes that the samples 
inhabit (see Table [3j> - We discuss the mean motions of M 
dwarfs as traced by the solar reflex motion, the velocity 
dispersions exhibited by active stars compared to their 
inactive counterparts and the change in velocity disper- 
sions seen in samples at different galactic heights. 

4.2.1. Solar Peculiar Motion 

The statistical parallax method measures the reflex so- 
lar motion with respect to the mean velocity of the stellar 
subsample being analyzed. If subsamples possess differ- 
ent mean velocities, this will be manifested as a change 
in the solar reflex motion. For reference, we use a co- 
ordinate system with U increasing towards the Galactic 
center, V increasing in the direction o f solar motion, and 
W in creasing vertically upward (as in IDehnen fc Binnevl 
1998). This system is right-handed, with the angular 
momentum vector of the solar orbital motion pointing 
towards the South Galactic Pole. 

In the U and W directions, the solar reflex motion re- 
flects the Sun's peculiar motion, as these distributions 



should be centere d on the local standard of rest (LSR; 
iFuchs et al.l 120091 ) . The solar reflex motion in the V di- 
rection is a combination of the Sun's pe culiar motion and 
the a symmetric drift at the solar circle (|Stromberdll924 
[19251 ). This shifts the mean V velocity of a sample of typ- 
ical disk-age M dwarfs, giving a V reflex motion larger 
than one measured from a population of young stars. 
Our results are shown in Figure [4j and we compare our 
peculiar solar velocities (i.e. the negative of the reflex 
motion) to previous work in Table [5] Our values were 
computed by taking the weighted mean across all spec- 
tral types with the uncertainty given by the standard de- 
viation. While the U and W velocities remain relatively 
constant with spectral type, the V reflex motion is sig- 
nificantly larger at early spectral types (M0-M2). These 
stars are observed at larger distances (Tabled, and also 
have larger velocity dispersions (as seen in Figure [5]). 
The same behavior is observed for the bluest r — z colors, 
which contain the same early-type M stars. For r — z > 2, 
the V velocity of the Sun remains relatively constant at 
20 km s _1 . This value compares favor ably to the V ve- 
locitie s derived by lHawlev et al.l ()1996l ) and IFuchs et al.l 
(2009), both measured with M dwarfs, but is larger than 
the V ~ 5 km s _1 value reported by IDehnen fc Binnevl 
( 1998) from measurements of the Sun's motion through a 
very young population, due to the age range of the stars 
that populate our subsamples. At later types, where the 
separation in age is most extreme between the active and 
inactive populations (see lWest et al.|[2008l and discussion 
in Sec 4.2.2 below), the active stars in the sample show 
decreasing V velocities, while the inactive stars have in- 
creasing V velocities, which supports the age hypothesis. 

As mentioned above, the measured solar reflex veloc- 
ity is the relative motion between the Sun and the mean 
velocity of a subsample of M dwarffl Since the Sun's 
velocity is not changing, a measured change in the reflex 
velocity indicates a difference in the mean velocity of a 
particular subsample with respect to the LSR. While the 
U and W velocity dispersions are expected to be constant 
with color (or spectral type), we find that there is sig- 
nificant structure in both distributions. The W velocity 
distribution in Figure [4] displays a non-monotonic behav- 
ior, peaking near r — z ~ 2.3 (type M4). This structure 
indicates the mean vertical motion of the particular M 
dwarf subsample is varying, with both bluer and redder 
M dwarfs having a smaller mean velocity. The W ve- 
locities of active and inactive stars are not significantly 
different. Meanwhile, the U velocity distribution shows 
a different behavior, remaining relatively constant from 
M0-M5 (r — z ~ 2.5), and then exhibiting a rise toward 
later types. The active stars appear to show this rise in 
U velocity at a bluer color (r — z ~ 2.0), but then are 
joined by the inactive stars at later types. This may be 
due to the changes seen in the other (V, W) velocities 
and/or to the change in the absolute magnitude of the 
active stars in bins where they are a significant fraction of 
the total distribution. Modeling of the populations to ex- 
plore changing the various parameters (activity fraction, 
absolute magnitude, input velocities in each direction) is 
necessary in order to disentangle these related effects. 

8 We note that the solar motion has also been estimated by 
observing th e direction a nd velocity of interstellar He entering the 
heliosphere i Witte! 120041) . 
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Fig . 2. — Left Panel: M r vs. spectral type. The open circles are nearby stars with accurate trigonometric parallaxes from B ochanskil 
(2008). Previous spectroscopic parallax relations are overplottcd and described in the legend. The results from this study (black, red and 
blue lines and black filled circles) are plotted for comparison. Note the large dispersion in M r for a given spectral type, especially at M4 
and M5. We strongly suggest that spectral type should not be used to estimate absolute magnitude in the SDSS r-band. Right Panel: M r 
vs. r — z, with same symbols and lines as in the spectral type panel. Previous studies are shown and described in the legend. There is a 
much smaller dispersion in M r at a given r — z color compared to the spectral type relations. 



10 



12 



14 



16 



i i i 



£ = 0.88 (Inactive) 
£ = 0.88 (Active) 
f= 1.1 (Inactive) 
, 5 = 1.1 (Active) 
All Stars 




_i i i i i_ 



10 



12 



14 



16 



| i i i i i i i i i | i i i i i i i i i | i 

£ = 0.88 (Inactive) 
£ = 0.88 (Active) 
£=1.1 (Inactive) - 
£=1.1 (Active) - 
All Stars - 




I 



M0 M2 M4 M6 
Spectral Type 



2 

r-z 



Fig. 3. — Left Panel: M r vs. spectral type for two values of a metallicity proxy. Active stars (dashed and dotted lines) and inactive 
stars (dot— dash lines) at the same metallicity are plotted. The total sample is overplotted in the solid black line. At a given metallicity, 
active stars are brighter than their inactive counterparts, while low-metallicity stars (f = 0.88) are dimmer than higher metallicity M 
dwarfs (£ = 1.1) at the same spectral type and activity state. Right Panel: The M r vs. r — z results show the same features as in the 
spectral type panel. 
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TABLE 6 
Solar Peculiar Motion 

Study £/© a V Q a ' b Vt© a 

Hawl ev et al. (1996) 9.1 ± 2 23.3 ±2 7.6 ± 2 

Dehnen fc B7nnevTl998') 10.00 ± 0.36 5.25 ± 0.62 7.17 ± 0.38 

Dehnen fc Binnev (1998) c 10 ±1 22 ±1 8 ± 2 

Fuchs et al. f2009) 9 ± 1 20 ±2 7 ± 1 

Aumer fc Binnev (2009) 9.96 ± 0.33 5.25 ± 0.54 7.07 ± 0.34 

Bond et al. (2010) 10 d 20 6.5 ± 0.4 

Schonrich et aTT^OlO) 11.1 ± 1 12.24 ± 2 7.25 ±0.5 

This Study c 8±2 24 ±3 7±2 



a UVW velocities are reported in km s , with U increasing 
radially inward, V increasing in the direction of the Solar 
orbit, and W increasing above the Galactic plane. 
b Measured relative to nearby M dwarfs, except for the 
Dehnen & Binnev (1998) and A umer fc Binnevl ((20091 ) stud- 
ies. 

c Estimated from reddest bin in Figure 3 of 
IDehnen fc Binnevl fl998l ). 
d Assumed. 

Error-weighted mean over spectral types M0-M7. 



4.2.2. Velocity Dispersions 

The active M dwarfs are known to form a kinemati- 
cally colder populat i on than their inacti ve counterparts 
(e.g., iWielenl fl977t lHawlev et all I1996D . as evidenced 
by smaller velocity dispersions. This is usually in- 
terpreted as an age effect, where younger M dwarfs 
have had less time to interact with giant molecular 
clouds and spiral density waves. The velocity disper- 
sion is predicted to vary as the square-roo t of age (e.g., 
iFuchs et al.ll200lHHanninen fc Flvnnll2002D . Young stars 



also harbor strong magnetic fields that heat their chro- 
mospheres, allowing actiy ity to be used as a proxy for age 
(| Wilson fc Woollevl I197C0 . Activity in low-mass stars 
has been shown to depend on both spectral type and 
age; inactive, late-type M dwarfs a re usually much o lder 
than inactive early-type M dwarfs (jWest et al.ll2Q08l ). In 
Table [71 we list the UVW velocity dispersions for the ac- 
tive and inactive subsamples by spectral type and color. 
Figure [5] shows the results. As expected, the dispersions 
of active stars separate from their inactive counterparts 
at the same spectral type (color) for mid-late M dwarfs. 
In the first bin where there are enough active stars to 
carry out the statistical parallax analysis (at spectral 
type M2), the active stars have larger reported disper- 
sions than the inactive ones. The velocity ellipsoid in 
this case shows a strong vertex deviation, indicating a 
young population, which has not undergone significant 
dynamical heating (Bi nnev fc M crrifield 1 !)()>). and the 
velocity dispersions should not be interpreted as mea- 
surements along the UVW axes. The active M2 stars 
being members of a young population is consistent with 
age-activity relations th at predict short ac tivity lifetimes 
for early-type M stars. (jWest et al.ll2008[ ). 



TABLE 7 
Velocity Dispersions of M dwarfs 



Spectral Type 




(km s- 1 ) 




(Tv- 


(km s 1 ) 




aw 


(km s" 1 ) 






Active 


Inactive 


All 


Active 


Inactive 


All 


Active 


Inactive 


All 


M0 




48 


47 




67 


67 




29 


29 


Ml 




45 


46 




58 


59 




29 


29 


M2 


59 


44 


45 


44 


66 


48 


34 


28 


29 


M3 


43 


47 


47 


29 


51 


42 


29 


29 


29 


M4 


36 


16 


44 


41 


11 


37 


24 


26 


26 


M5 


31 


19 


42 


23 


39 


39 


18 


25 


22 


M6 


27 


39 


33 


17 


28 


27 


13 


23 


18 


M7 
r — z 


27 


13 


34 


17 


29 


28 


12 


23 


17 


0.75 




49 


48 




66 


65 




26 


26 


1.05 




47 


47 




63 


63 




30 


31 


1.35 




42 


44 




56 


57 




27 


28 


1.65 


47 


45 


46 


39 


42 


44 


33 


27 


27 


1.95 


41 


50 


48 


60 


11 


44 


28 


28 


29 


2.25 


34 


45 


43 


24 


36 


34 


19 


27 


26 


2.55 


33 


47 


41 


26 


29 


28 


18 


22 


21 


2.85 


26 


12 


36 


17 


29 


25 


12 


22 


18 


3.15 


27 


39 


32 


17 


26 


23 


12 


22 


16 


3.45 


27 




32 


13 




20 


11 




17 



Note. — The typical uncertainty for the dispersions is 5-10 km s 
Equation 



Another notable feature of Figure [5] is the decline in 
velocity dispersion at later types. Since the type (color) 
bins sample different spatial volumes (see Table [3]) , this 
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Fig. 4. — Solar peculiar motion as a function of spectral type (left panel) and r — Z color (right panel). The V velocities exhibit a decline 
to later type (or redder) stars, likely due to an age effect, as shown by the separation between the active and inactive populations. The 
W velocities exhibit a peaked behavior, while the U velocities show a rise toward later types, which occurs at a bluer color in the active 
population. 
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Fig. 5. — UVW velocity dispersions vs. spectral type (left panel) and color (right panel) for active (red triangles), inactive (blue squares) 
and all stars (filled black circles). Active stars possess smaller dispersions at mid-late spectral types (red colors). The M2 bin shows a 
strong vertex deviation so the velocity dispersions should not be interpreted as measurements along the UVW axes. 



is likely showing that stars that are located further above 
the plane of the Galaxy (in the early-type bins) are there 
because they are older and have undergone more dynam- 
ical interactions, giving correspondingly larger velocity 
dispersions. To investigate this kinematic structure in 
the Milky Way, we use the IPA parameter described in 
Equation [3] as a proxy for height. 



The mean height above the Plane, Z, can be approxi- 
mated as: 

Zwdsinfe (4) 

where b is the mean Galactic latitude for the subsample 
(~ 50° for SDSS as a whole). This approximation ne- 
glects e ffects due to the Sun' s location above the p lane 
(15 pc; lCohenlll995l ). but see (|Bochanski et al.ll2010l) for 
a more precise expression. Given the definition of IPA in 
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Equation [31 Z and IPA are thus related by: 

Z fn IPA lcK 1- ^). (5) 

For reference, we tabulate various heights for a given IPA 
and M r in Table |H To increase the number of stars in 
each IPA bin, we used stars at the same absolute distance 
from the Plane in both the northern and southern Galac- 
tic hemispheres, thus implicitly assuming the Galaxy is 
symmetric about its mid-plane. 

TABLE 8 

Heights for various IPAs and Absolute Magnitudes 81 



IPA 












8 


10 


12 


14 


100 


25 


10 


3 


1 


1000 


251 


100 


39 


15 


3000 


753 


300 


119 


17 


5000 


1255 


500 


199 


79 


6000 


1507 


600 


238 


95 



a All heights reported in parsecs. 

Pre v ious k in ematic inve s tigati on s (iBochanski et alj 
2001 I2007at IWest et all [20081 iFuchsetalJ 1200a 
Schmidt et al.l 120101 ) have demonstrated an increase in 
velocity dispersion at larger Galactic heights, rising from 
((7(7, cry, <Jw) ~ (30 km s _1 , 20 km s~ x , 20 km s _1 ) in 
the local solar neighborhood, to (cry, oy , <Jw) ~ (50 km 
s" 1 , 40 km s~ 1 , 40 km s _1 ) at 1 kpc. In Figure [51 we 
plot the vertical W velocity dispersion as a function of 
IPA (left panels) and mean height above the Plane (right 
panels), for spectral type (upper row) and color (lower 
row). Recall that the advantage of the IPA is that no 
absolute magnitude is assumed, while the heights are 
calculated using the absolute magnitude from the sta- 
tistical parallax solution for a given subsample. The in- 
crease in velocity dispersion is seen in both height and 
IPA. Furthermore, there is no strong dependence on spec- 
tral type or color, as expected since the dynamical in- 
teractions should be controlled by the Galactic gravi- 
tational potential, with the stars (of any mass) acting 
as collisionless test particles. These results bridge the 
gap b etween the kinematic investigations of [Bond et al.l 
(2010), which examined stars beyond our distances lim- 
its, an d those examining nearby stars interio r to our 
study ([Nordstrom et al.|[200l iReid et al.ll2002ft . 

The velocity dispersions for M dwarfs listed in Table 
[7] should be interpreted in light of these IPA/height re- 
sults. Earlier spectral types (i.e., M0 and Ml) sample a 
larger range of distances and have a larger mean distance 
than later types (see Table [3j . Figure [5] shows that both 
early-type and late-type M dwarfs have similar velocity 
dispersions when sampled at the same Galactic height. 

5. CONCLUSIONS 

We present a statistical parallax analysis of the most 
recent ca talog of M dwarfs identified with SDSS spec- 
troscopy ([Paper J ). Our sample was subdivided on many 
criteria, to explore both the intrinsic changes in low-mass 
dwarfs due to metallicity and magnetic activity, and to 
investigate their kinematics. We have demonstrated that 



r — z color is a much better proxy for absolute magni- 
tude than spectral type, and suggest that photometric 
parallaxes are the preferred method to determine abso- 
lute magnitudes and distances for M dwarfs. 

Some interesting trends with M r were revealed in our 
analysis. First, magnetically active M dwarfs were shown 
to be intrinsically brighter at the same spectral type 
or color than their inactive counterparts. Eclipsing bi- 
nary studies have dem onstrated that magnetic activity 
inflates a star's rad i us (IRibasll2006t IL6pez-Moralesll2007l; 
IMorales et al.ll2008ll2010[ ). but the effect on effective tem- 
perature and luminosity is less constrained. Our results 
suggest that activity may increase the radius and lumi- 
nosity of active low-mass stars. 

Metallicity, which influences the luminosity of a star, 
was also ex plored. We divide d the sample using the Q 
parameter ([Lepine et al.1 l20"07f) as a proxy for metallic- 
ity. Our results are similar to those observed for higher 
mass stars: low-metallicity M dwarfs are dimmer at a 
given spectral type (or color) than high metallicity stars. 
We have quantified this effect for M dwarfs in the SDSS 
photometric system. 

To isolate the effects of metallicity and activity, we 
separated active and inactive stars for the same ( as a 
function of color and spectral type. Activity still bright- 
ened stars at the same £, however the effect diminished 
at smaller £ (lower metallicity). 

The statistical parallax analysis also allowed us to in- 
vestigate the reflex solar motion and velocity dispersions 
for each subsample. The more distant, early-type stars, 
which are presumably older, have the largest reflex so- 
lar motion, particularly in the V direction, which we at- 
tribute to the increased asymmetric drift. The active and 
inactive stars exhibit expected behavior with the active 
star populations having smaller mean motions relative 
to the Sun. The inactive, late-type M dwarfs, which 
we identify as an older p opulation due to their lack of 
activity (ages > 4 Gyrs, IWest et all 12008ft have veloc- 
ity dispersions similar to early-type M dwarfs, which we 
identify as old due to their greater vertical distance from 
the Galactic plane. Thus, the activity and dynamical 
heating age indicators give consistent results. When the 
velocity dispersions are analyzed as a function of vertical 
height or independent position altitude (IPA), all stars 
exhibit increasing dispersion at increasing height above 
the Plane. 

As astronomy enters a new e ra of large photom etric 
survey s, such as PanST ARRS ( Kaiser et all 12002ft and 
LSST (|Ivezic et al.ll2008ft . it will be vital to develop tech- 
niques for estimating the activity and metallicity of low- 
mass stars from photometry alone. There have already 
been efforts to characterize metallicity usi ng the color- 
color distributions of M dwarfs in SDSS (Lepine 2009; 
iPaper J ), but those data must be calibrated with spec- 
troscopic observations. The work presented in this paper 
highlights the need to determine these fundamental pa- 
rameters, since they affect the estimated distance to each 
star. Finally, the importance of the dual spectroscopic 
and photometric nature of SDSS cannot be overstated. 
The large spectroscopic samples of M dwarfs it has ac- 
quired have enabled many novel investigations, including 
this one. Significant spectroscopic followup of the next 
generation of surveys should be a high priority. 
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Fig. 6. — W velocity dispersions vs. IPA (left panels) and height above the Plane (right panels). The velocity dispersions increase with 
IPA/height but do not depend on color/spectral type. 
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